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Abstract 

A class of singular integral operators, encompassing two physically relevant cases 
arising in perturbative QCD and in classical fluid dynamics, is presented and an- 
alyzed. It is shown that three special values of the parameters allow for an exact 
eigenfunction expansion; these can be associated to Riemannian symmetric spaces 
of rank one with positive, negative or vanishing curvature. For all other cases an 
accurate semiclassical approximation is derived, based on the identification of the 
operators with a peculiar Schroedinger-like operator. 
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1. Introduction 

It has recently been realized that a special kind of singular integral equa- 
tion arising in the study of jet production (e + e~ — » qq + anything) |2j 
bears a striking similarity to another equation 4 introduced fourty years ago 
by E. Tuck jS] in the context of laminar flows around slender bodies. In this 
note we describe a general two-parameter family of integral operators which 
reduce to Tuck's and Marchesini-Mueller's (hereafter mm) case for special 
values of the parameters: 

(v mm \ - f 1 ^0) - <Ky) i , 

{K a /3P){x)= / — ; ; — dy + 

J-i \x-y\ 

+ ((1 - a) log(l + x) + (1 - (3) log(l - x)) 4>{x) . 

We shall show that the spectral problem can be solved exactly in three cases, 
connected to the three distinct symmetric spaces of rank one (with curva- 
ture 1, 0, -1) and they correspond to (a, (5) = (1, 1), (0, 0), (0, 1), respectively. 
The first case [3] has discrete spectrum and it is (unitarily equivalent to) 
a function of the Laplacian on the sphere restricted to the axially symmet- 
ric sector. The second case has continuous spectrum and it is (unitarily 
equivalent to) a function of the Laplacian on the real line. The third case is 
equivalent to Marchesini-Mueller's operator and it is (unitarily equivalent 
to) a function of the radial Laplacian on the hyperbolic plane. Tuck's case 
corresponds to T = \K22 and no exact solution is presently known. To it 
we can nonetheless apply a semiclassical approximation (WKB) which will 
be derived in general for any positive value of the parameters; the spectrum 
is purely discrete in this case and it is approximated by 

(1) **-»h*+»>-*(ra$)- 

+ (l-±(a + /?)) log2 + 7 £ 

(je is Euler's constant), which in particular gives Tuck's eigenvalues to a 
very good accuracy (see Table 1). 

This note is organized as follows. In Section |2] we introduce the special 
problem related to jet-physics and we show that it is unitarily equivalent 
to Kqi. In Section El we identify a second order differential operator C 
commuting with Kq i and determine its eigenfunction expansion. Moreover, 
a first order differential operator £ is shown to commute with the operator 
Kqq, and also in this case we can obtain the spectral representation which 
is used in Section 0] to introduce another representation of the operators, 
equivalent to Schroedinger's equation with a kinetic energy given by g(p), 
where p is the momentum operator — id/du and g is essentially Lipatov's 
function. In this representation it is easy to derive qualitative properties 



We are indebted to R. A. Askey for pointing out to us the reference to Tuck's paper. 
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of the operator K and to set up the semiclassical approximation. We also 
derive the boundary behaviour of eigenfunctions in the general case. In 
Appendix A we show that the operator C commuting with Kq\ is indeed 
equivalent to the Laplace operator on the hyperbolic plane, a fact which 
gives us valuable information on the eigenfunction expansion (completeness, 
spectral measure). In Appendix B we give an easy proof of Tuck's result 
about K\ i which is essential to the developments of Sec. El 

2. Marchesini-Mueller's equation 

Marchesini and Mueller [Q, introduced an equation for the multiplicity 
of quark-antiquark pairs in electron-positron collisions. As a function of 
energy, the multiplicity density satisfies an evolution equation given by 



du(r, £) _ f 1 drj 



dr J 1 - 77 



1 — 77 



77 

where t is the logarithm of the energy in the center of mass and £ = 
^(1 — cos#), with 9 the angle between the two jets emerging from the 
electron-positron collision. Knowing the multiplicity at low energy, its en- 
ergy dependence can be calculated at all higher energies by QCD perturba- 
tion theory and the result, in a special regime, gives the integral equation 
above. For details see It is formally very similar to the so-called BFKL 
equation The unknown n(r, £) is defined for £ £ (0,1), and it vanishes 
at £ = to ensure convergence. The initial value problem is solved if we can 
find the spectral decomposition of the operator in r.h.s. 

It's a matter of simple algebra to show that actually the equation can be 
recast into the form 

-7^7 = ~Kq 1 (f) 

by performing the following transformation 

u(t,0 =e- lo s 2T £0(T,2£-l) 

The operators K\ \ and Kq 1 play a central role in the following. Hence 
we introduce a special notation for them: 

H = ^Kn 

M = ^ i-log2 = 2W + logi(l + x). 

It is known [S] that TC is diagonal in the basis of Legendre polynomials and 
its discrete eigenvalues are given by the the harmonic numbers 



HP n = t) n P n , t) r 




for n = 
for n > 



(a simple proof of this result can be found in Appendix B). One could study 
the general spectral problem for (a, /3) close to (1,1), e.g. by perturbation 
theory. However for arbitrary values of the parameters a different approach 
is needed. It has been shown in Ref.[2] that the operator M has actually a 
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continuous spectrum, and the eigenfunctions can be identified with hyper- 
geometric functions. The result is obtained by an expansion starting from a 
combination of phase-shifted plane waves. The expansion can be pushed to 
all orders, and the resulting series is convergent to a hypergeometric function 
which can be identified with Legendre functions. The spectral decomposi- 
tion of M is reduced to the classical Mehler-Fock transform. We shall come 
back to these facts in Appendix. Here we want to show how this exact result 
can be derived without any approximate procedure, by looking for a local 
(differential) operator commuting with M. This will give an alternate more 
rigorous proof of the solution. 



3. COMMUTATIVITY WITH DIFFERENTIAL OPERATORS AND EXACT 

SOLUTION 

3.1. The mm operator. The easiest way to solve mm equation is to find a 
differential operator C that commutes with M: 

[£, M] = [C, 2H] + [C, log(l + x)]=0 

It is convenient to look for the operator C such that C as well as [£, log(l+x)] 
act in a simple way on P n (x); taking into account the known properties of 
Legendre polynomials 

£ P n = [{l-x 2 )d 2 x -2xd x ]P n (x) = -n(n+l)P n (x) 
n + 1 n 

n = 2^TT " +1 + 2^TT n_1 
-(i- x 2 )d x p n = !^±11 ( P n+1 _ P n _ l} , 

the problem will reduce to a purely algebraic one. Since the mm equation 
has a singularity at x = — 1 we will search C in the form: 

C = (1 + x) C + a(l - x 2 ) d x + 6(1 + x) 

By construction, this operator acts in a simple way on P n {x). Namely, 

£ P n {x) = A n P n +l + B n P n + C n -i P n -\ 

where 5 

j(2n+l)A n = -n(n + l) 2 - an(n + 1) + b (n + 1) 
\ (2n+ l)C n „i = -n 2 (n+ 1) + an(n + 1) + bn 



(2) 



(and B n does not enter in what follows). The action of the commutator 
[jC,TC] on P n is then given by 

(3) [£, H] P n = hnCP n - H£P n = -C n _iP n _i - -^—A n P n+1 

n n + 1 



5 Our convention for the indices is the natural one if we think to the r.h.s. as the action 
on the left by a tridiagonal matrix with vectors [C, B, A] along the diagonal. 
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The commutator C = [C, log(l + x)] can also be easily calculated: 

(4) C = 2(1- x 2 ) d x - 2x + {a - - x) 

Notice that there are no diagonal terms coming from Eq. (|5|>. hence the 
only diagonal contribution to the commutator comes from the last term in 
C, which immediately implies a = 1. Using the properties of P n , we can 
write: 

CP n = R n P n +i + S n -i P n -\ 

with 

{n + lf _ o n 2 



2n + l 2n + l 

Now we can check the commutativity of £ and A4: 

2 2 

+ Rn = 0; — C n _i + 5 n _i = . 



n + 1 n 
From the first equation we have 

(5) A --^±1r - (n+1)3 - 

i5j An ~ 2 " 2n + l ' 

on the other hand we have from Eq. © 

(n + l)(n 2 + 2n - b) 
n ~ 2n+ 1 

which fixes 6 = — 1. The equation for C n -\ is automatically satisfied. Notice 
that for K a ± with a > one would find a coefficient different from one in 
front of A n in Eq. (JSJ), hence no solution. A4 is therefore the only operator in 
the family K a g which allows a commuting differential operator of the form 
C The other two cases, alluded to in the Introduction, are connected to a 
different choice of £, and will be discussed later. 

Now we can find the eigenfunctions for the mm equation. They satisfy 
the differential equation: 

£0 = \<f>. 

We should look for the solutions satisfying \<j)(x)\ — > |1 + x\~ 1 / 2 at x — > — 1 
and which are finite at x — ► 1. It is convenient to parametrize A = —1/2 — 
2k 2 . Then the equation for (j> can be written as: 

1 -l- Ah 2 

[(1 _ x 2 )d 2 + (1 - 3*R - l + ^±_^ = 

With the substitution: 

(j) = (l + x ) a ^(x); a = -1/2 + ik 

we get 

[(1 - x 2 )d 2 x + {2a + 1 - (3 + 2a)x}d x (a + 1) 2 ]V>0) = . 
Setting x = 1 — 2y, we obtain: 

[y(l - y)8 2 + {1 - (3 + 2a)y}d y - (a + l) 2 ] V(y) = 



which is a hypergeometric equation with a = b = a + 1 = 1/2 + ik and c = 1. 
The solution which is finite at y = (x = 1) and at y = 1 (cc = — 1) has the 
form: 

i/j(y) = F(l + ik,l + ik,l,y); 

</>{k,x) = (l + x)-^ +ik F(\ + ik,\ + ik,l,^^) 

If we re-introduce the variable £ = (l + ic)/2 the result can be written in the 
form: 

(6) <p(k, o = c r 1/2+ik f{\ ik, 1,1-0 

To calculate the eigenvalue k of MM as a function of k (the dispersion rela- 
tion) we may use the fact that H annihilates the constant and it symmetric. 
It follows 

(7) K (k) [ </>(k,£)d£= [ log^(fc,0d£ 



This integral can be calculated and it is given by a combination of digamma 
functions known as Lipatov's function ^ 

= ^(1 + ik) + - ik) - 2 ip(l) 

= -4 log 2 + 14C(3) k 2 - 62C(5) A; 4 + 0(k 6 ) , 

(il>(z) = dlogr(z)/dz). The evolution equation for mm equation can now be 
solved by expanding u(r, £) on the continuous basis £ The spectral 

measure which defines the eigenfunction expansion can be taken by [3] (see 
also Appendix). 

3.2. The case a = f3 = 0. We note that another operator, namely: 2TC + 
log(l — x 2 ) = Kqq also commutes with a differential operator. In this case 
it is easy to prove that it commutes with a first order differential operator 6 



(9) £ = i[-(l - x 2 ) d x + x] = -n/l - x 2 d x y/T- x< 
The eigenfunctions in this case have a simple form: 

£4>(k,x) = kcj)(k,x) 

(10) <f>{k,x) = (1 + x) {ik ~ 1)/2 (1 - x)- {lk+1)/2 . 

The eigenvalue g(k) of Kqq = 27i + log(l — x 2 ) belonging to these eigen- 
functions can be calculated exactly in the same way as before (see Eq. (J7J)) 
and it has the form: 

g(k) = V>(|(1 + ik)) + ^(|(1 - ik)) - 2 V(l) + 21og2 

(which again can be reduced to Lipatov's function). This means that Kqq 
coincides with a function of the first order differential operator, namely 

(11) 2H + log(l - x 2 ) = g(£) 



6 It will be noticed that this operator is a multiple of the commutator C of Eq.@. 
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In the next section we will use this representation to derive a semiclassi- 
cal approximation and for the analysis of the asymptotic behaviour of the 
solutions (j)(x) near the boundary points x = dbl. 

4. The Schroedinger representation and semiclassical analysis 

4.1. The semiclassical spectrum. We can use the representation (fTTj) for 
the operator 27i + log(l — x 2 ) in terms of the first order differential operator 
to transform our integral equation in the form of a Schroedinger equation 
(with unusual kinetic term) which is convenient for the semiclassical analysis. 
Consider the equation: 



(12) K al3 cj ) = {2H + l(l-a) log(l+ x) + (1-/3) log(l- x)]} 4> = k<\> . 



This equation can be rewritten in terms of Schroedinger equation. Namely if 
we do the substitutions: x = tanhn and <p = cosh(it)\I/(ii) the last equation 
can be rewritten as: 

(13) g(-id u )^(u) - [alog(l + tanhu) + /?log(l - tanhn)]^ = 

In the free case a = (3 = we have the plane waves solutions corresponding 
to the functions (|1U|) after this substitution. In the case a = {3 = 1 we have 
'I'll = Pn (tanh u) I coshu; this means that, by identifying tanhu = cos??, 
and modulo a similarity transformation, this case is related to the Laplace 
operator on the two-dimensional sphere. Finally, the case a = 0, {3 = 1 
corresponds to MM equation, where also the exact solution is known. 

It is convenient to slightly modify the function g{z) and the eigenvalue n 
by adding a constant shift. Let us introduce G(z) and «' by: 



The operator G(p) takes on the role of the kinetic energy and is equivalent to 
the usual operator of non-relativistic quantum mechanics in the low energy 
limit. It follows that the operator G(p) + V(u) has a discrete spectrum for 
(3 ^ a > (by symmetry we can restrict to the sector (3 ^ a with no loss of 
generality), while it has continuum spectrum for (3 ^ a = 0. Let us notice 
that in this picture Tuck's operator T (a = (3 = 2) is qualitatively very 
similar to the "trivial" case {a = (3 = 1) which corresponds to operator 7i 

The "Schroedinger" representation can be used in different ways. In par- 
ticular, we can exploit this representation to compute the semiclassical ap- 
proximation to the eigenvalues and eigenfunctions. We note that near the 




We shall need the following asymptotic behaviour of G(p): 
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turning points, where the kinetic term is small we can apply the standard 
WKB approach. This leads to the following (Bohr-Sommerfeld) semiclassi- 
cal approximation to the eigenvalues: 
r b 

(17) / G-\k' -V)du = 7r(n + l/2) 

J a 

where G~ l is the inverse function of G and a, b are the turning points. It 
follows from Eq, (|16|) that the inverse function is approximately G _1 (£) = 
exp(£/2) + 0(exp(— £/2)). In the region where the semiclassical approxima- 
tion works we can neglect all asymptotic terms besides the first one. In the 
main approximation we can also put a = — oo, b = +00. Then equation (|17|) 
can be rewritten as: 

+00 

(18) J exp(i( K ' - V))du = vr(n+ 1/2) 

—00 

This integral can be easily calculated and we derive for the spectrum of the 
operator K a p\ 

4 a '« « 2 [log(vr(n + 1/2)) - log (B(a/2, (3/2)) + (1 - \ (a + /?)) log 2 + lE ] , 

where B is Euler's "beta" function i.e. we arrive at Eq. This equation 
gives the standard approximation for the harmonic numbers f) n (the eigen- 
values of \K\i) up to 0(l/n 2 ), while for the eigenvalues of Tuck's operator 
T = \-K22 it gives: 

(19) ±4 2 ' 2 ) « log(7r(n + 1/2)) - log 2 + 1E . 

The semiclassical formulas are contrasted to the numerical or exact eigen- 
values in Tab.l for operators T and H. respectively. 



n 


Ref. A 


WKB 




WKB 





0.2332 


0.3357 





-0.116 


1 


1.4437 


1.4343 


1 


0.9827 


2 


1.9409 


1.9451 


1.5 


1.4935 


3 


2.2833 


2.2816 


1.8333 


1.8300 


4 


2.5317 


2.5329 


2.0833 


2.0813 


5 


2.7342 


2.7335 


2.2833 


2.2820 


6 


2.9000 


2.9006 


2.4500 


2.4490 


7 


3.0440 


3.0437 


2.5929 


2.5921 


8 


3.1686 


3.1689 


2.7179 


2.7173 


9 


3.2803 


3.2801 


2.8290 


2.8285 



Table 1. The WKB spectrum of T and H. 
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The semiclassical eigenfunctions ^f^ sc '(u) can be written in the form: 
where A is a normalization factor and K f n — K n — 

2^e- The integral gives the 
incomplete Beta function which reduces to elementary trascendentals in the 
special cases cases a = f3 = 1 and a =(3 = 2 corresponding to the operators 
7i and T respectively. In the first case we have: 

(20) = Al sin[(2n + l) tan-He") + ,/4] 

V cosh u 

with = (±tt + 2^h) _1/2 - In the 

case corresponding to Tuck's operator 

T we obtain: 



(21) ^ c \u)=A 



sin[7r/2(n + l/2)(tanhn + 1) + vr/4] 



coshu 

with A2 = (1 + 2n+i )~ 1 ^ 2- ^ we rewr ite these functions in terms of original 
variables x = tanhu and functions 4>(x) = cosh(u)\I / (u) then the expression 
H2U|) gives the well known large n asymptotics of Legendre polynomials: 

sin[(n + 1/2)0 + 7T/4] 
n (cos6») r — = . 

V sin t) 

For Tuck's case (a = (3 = 2) the semiclassical wave functions have a very 
simple form: 

(22) <j) n (x) =,4sin[7r/2(n + l/2)(x + l) + 7r/4] 

and they give a rather accurate description (Fig.l) of the true eigenfunc- 
tions which can be easily computed numerically {i.e. by using the spectral 
representation for the operator 7i on the Legendre basis). 

4.2. Boundary behaviour. The "Schroedinger" representation (|15|) can 
be used to derive the asymptotic behaviour of the eigenfunctions 4>{x) at 
the singular points x — * ±1. Namely we show that these asymptotics have 
the form: 

(23) <j) a p ~ I log(l + x)\ da I log(l -x)\ d P 
where 

(24) d a = --l; d = ^-l. 

a p 

We note that the asymptotic behaviour of the function 4>(x) at x — > ±1 
can be derived from the asymptotics of the function ^(u) at u — * ±00. 
Consider, for example, this asymptotics at u — > +00 (x — * 1). In this limit 
the potential term V has the form: 

V = 2(3u - 2a log 2 + 0{e~ 2u ) 
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We can neglect all terms in this expansion besides the linear one. In this 
way the problem is reduced to the calculation of the asymptotics at u —* oo 
of the solution of Eg, 1)15(1 with a linear potential. For this purpose it is 
convenient to rewrite this equation in the Fourier representation. It follows 
from the explicit form of function G(p) that the Fourier transform *$>(p) of 
ty(u) satisfies the first order differential equation: 

-2f3id p ^(p) + [i/>(l/2 + ip/2) + ^(1/2 - ip/2) + log4]*(p) = 
The solution of this equation has the form: 

(25) *(p) = (jr* r(|[i + ip)) ) ex P(^W2/3) 

The asymptotics of the function 

(26) ^(u) = J exp(-ipu)fy(p)dp 

is determined by the nearest singularity of the function ^>(p) in the lower half 
plane at the point p = —i. For non-integer 1/(3 this singularity is a branching 
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point. The standard estimation of the corresponding contribution gives: 

-^u d ^e~ u (l + 0(l/u)); u^oo 

Taking into account that x = tanhn and (f>(x) = cosh(ti)^(n) as well as the 
a, (3 symmetry of the equation we arrive to Eas. ()23l 124(1 . It follows from 
Ea.(|25|) that this asymptotics takes place in rather narrow region u ^> k' /2(3 
or | log(l — x)\ 3> k' / (3 (for x — > 1). 

We note that for f3= 1 the integral Q26|) can be calculated explicitly and 
we can derive the exact wave functions in the potential V = 2u. They have 
the form: 

=yJ (2y) 

where y = exp(— u + k'/2) and Jq{z) is the Bessel function. 
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Appendix A: The connection to hyperbolic plane 

The eigenfunctions 4>(k, x) given by Eq. © are related to Legendre func- 
tions with complex index. To see this, let us revert to the original form of 
MM equation for which the eigenfunctions are given by 

u(k, x) = x 1/2+tk F{\ + ik,\ + ik, 1, 1 - x) . 

By using well-known properties of Jacobi functions (see 0) we find that 

u(k,x) can be identified with the Jacobi function <^2fc°^ (*) W1 th x = 1/ cosh 2 (t) 
Hence it follows that u(k,x) = P_i/2+ik{2/ x — 1). Now, it is well known 
that Legendre functions of this kind appear as spherical functions on the 
hyperbolic plane , i.e. they are eigenfunctions of the radial part of the 
Laplace operator in the case of the two-dimensional homogeneous Riemann- 
ian space with constant negative curvature. In Gaussian coordinates, ds 2 = 
dr 2 + sinh 2 r dip 2 , the radial part of the Laplacian A r is given by 

A r = (d/dr) 2 + coth r d/dr . 

and it is immediate to check that (A r + (1/4 + A; 2 )) P_ 1 ^ 2 +ifc(coshr) = 0. 
It is then natural to conclude that there must exist a map x — > r and a 



7 see e.g. 00. 
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suitable similarity transformation which connects the operator C of Sec. 01 
to A r . Prom the expression of u(k,x) in terms of Legendre functions, the 
map is given by coshr = 2/x — l = 4/(z + 1) — 1, where z is the variable 
which enters the definition of C The similarity transformation is simply 

\C = (1+coshr) A r (l+coshr) -1 . 

Having established this connection, the explicit eigenfunction expansion 
comes for free, in terms of Mehler-Fock transform 

' oo 

u(x) = f P_ 1/2+lk (| - 1) c(k) dk, (0 <x <1) 

< 

c(k) = k tanh-Kk f u ( ^ ) P_ 1 / 2+i fe(t) dt 
\ l v y 

which can be used to solve the evolution in r for the mm equation 2 . 



Appendix B 

We give a simple proof of an old result due to E. Tuck 

Theorem. 7^11 has the Legendre polynomials P n (x) as eigenvectors with 
eigenvalues the harmonic sums h n . 

Proof. By computing K\\p n with p n (x) = x n we find 

/l X n — y n 
dy-j r 
-1 F-y 



/a; /• 1 \ n 

- / ) dl/ E ^ _1 ^~ fc 
- 1 / fc=l 

2n,..r" E' ! ""^ 



fc=i 

hence i^i 1 leaves each subspace P n of polynomials of degree n invariant for 
any n. It's matrix representation is upper triangular and its eigenvalues 
are found on the diagonal by inspection. Since K\\ is symmetric with re- 
spect to the inner product {p\,P2) = j\ dxp\(x)p2(x) its eigenvectors are 
orthogonal, hence they are the Legendre polynomials. □ 
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